Probing the Quantumness of Channels with Mixed States 
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We present an alternative approach to the derivation of benchmarks for quantum channels, such 
as memory or teleportation channels. Using the concept of effective entanglement and the ver- 
ification thereof, a testing procedure is derived which demands very few experimental resources. 
The procedure is generalized by allowing for mixed test states. By constructing optimized measure 
and re-prepare channels, the benchmarks are found to be very tight in the considered experimental 
regimes. 
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I. INTRODUCTION 

The field of quantum information science has brought 
forward a range of interesting information processing pro- 
tocols, many of which have already been demonstrated 
experimentally in a multitude of physical implementa- 
tions. In quantum communication, the most promi- 
nent examples are (arguably) quantum key distribution 
S]j teleportation [1, 0, 0, and quantum memory 
3, l8[ experiments, all of which were implemented in dis- 
crete (qubit) as well as continuous- variable settings. 

All of these protocols are based on the fact that in- 
formation is transmitted in the form of quantum states, 
either by directly sending non-orthogonal states or by us- 
ing pre-shared entanglement. This stands in contrast to 
measure and re-prepare channels, which convert incoming 
states to classical data by means of measurements and re- 
prepare new quantum states accordingly. This conversion 
to classical data destroys any advantage which quantum 
protocols have over classical ones. 

This observation leads to a natural benchmark cri- 
terion: An experimental implementation of a quantum 
channel is successful only if it can outperform any mea- 
sure and re-prepare channel. How to quantify the perfor- 
mance of a channel is not trivial. For memory or telepor- 
tation channels, this is typically done with the average 
fidelity 0. For a given ensemble of test states {Pi,p\ n }, 
and corresponding output states p° ut , the average fidelity 
is defined as F = T li Pi F (Pi l : P° ut )- 

For continuous-variable protocols, a maximization of 
the average fidelity for measure and re-prepare strate- 
gies was first performed by Hammerer et al. for a Gaus- 
sian distribution of coherent test states Narniki et 
al. recently derived an extension for channels of non- 
unity gain [ll|. Due to recent interest in the stor- 
age of squeezed light d, [Hj], fidelity-based benchmarks 
were derived for different ensembles of squeezed states 
[l3l [II . Hrl [Tfij . Clearly, finding the optimal measure 
and re-prepare channel is a challenging task even for in- 
put ensembles, which are very small or highly symmetric. 
Moreover, each benchmark is tied to an input ensemble 
in the sense that a change in the input ensemble requires 
a new optimization of the average fidelity. In addition, 



the fidelity can be hard to measure in an experiment, 
especially in the continuous-variable setting. 

In this paper, we use a different technique to derive 
benchmarks which are particularly simple to implement 
experimentally. Our analysis mainly relies on two ob- 
servations: Firstly, any measure and re-prepare scheme 
acts as an entanglement-breaking channel [17| . There- 
fore, it suffices to demonstrate that a memory or tele- 
portation setup can preserve entanglement. Secondly, 
every source of non-orthogonal states admits a theoret- 
ical entanglement-based description, which is typically 
referred to as source replacement or as effective entangle- 
ment [13, [l9|, l20( . This means that we can devise bench- 
marks based on the verification of entanglement without 
the need to generate actual entangled states in the labo- 
ratory. 

A general benchmark criterion should accommodate 
mixed test states, to cover cases where the generation of 
pure test states is experimentally infeasible. For exam- 
ple, initially pure test states may be subject to a noisy 
environment before they reach the quantum channel to 
be tested. Another example of this is the generation of 
squeezed light, where losses in the preparation process 
lead to mixing. We therefore extend the technique of 
source replacement to sources which emit mixed states. 

For the verification of effective entanglement, there ex- 
ists a large number of methods, each with their partic- 
ular advantages [2l[. Here, we will use the so-called ex- 
pectation value matrix (EVM) criterion [22|, [H| , which 
is tailored for situations where state tomography is not 
available. We show how to extend this criterion to incor- 
porate sources of mixed states. The methods are exem- 
plified with the help of three different physical settings: 
a qubit protocol, a source of displaced thermal states and 
a source of squeezed thermal states. 

This paper is structured as follows. Section |TT] gives 
a review of the concept of effective entanglement for 
sources of pure states. Furthermore, the expectation 
value matrix method is reviewed and applied to a source 
of squeezed vacuum states, which may be used to bench- 
mark quantum memory experiments. Section HTTl presents 
a formalism for source replacement and entanglement 
verification for mixed test states. Using the results from 
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the pure test states, this formalism is developed in steps, 
which approximate the actual problem increasingly well. 
In Sec. IIVI the effects of entanglement-breaking channels 
are considered. By optimizing the possible measurement 
outcomes induced by such channels, the whole domain 
compatible with separable states can be identified. This 
is a tool for assessing the strength of the developed en- 
tanglement criteria. Finally, Sec. fVl contains a conclusion 
and a discussion of the results. 



II. PURE TEST STATES 

A test for quantum channels using pure test states was 
proposed in Ref. [12] • It relies on an entanglement-based 
source description and the subsequent verification of en- 
tanglement. In this section, we start with an outline of 
the proposed test procedure and a brief review of the 
entanglement verification method. We then investigate 
the application of the method to quantum memories for 
squeezed light. 

The general setup is as follows: To probe the quan- 
tum channel, we employ a source of (non-orthogonal) 
test states \ip" 1 ), which are chosen at random from a 
fixed set {|V'l"}}t=o °f cardinality N. Each state is sent 
with a predetermined probability In quantum key 
distribution, this ensemble of test states is determined 
by the protocol, since it must be possible to exclude 
intercept-resend attacks using the same classical data 
that will ultimately be used to generate secret key. For 
quantum memories, the input ensemble can be chosen 
freely, in accordance with experimental feasibility. The 
key to the test procedure is the fact that every source of 
non-orthogonal test states admits an entanglement-based 
description. In this thought setup, an entangled source 
state \ip src ) = Yli=o y/pl\i) aW 1 ™) b 1S generated, and we 
label the two subsystems by A (Alice) and B (Bob). The 
correct test state ensemble can effectively be prepared by 
measuring system A in the basis {| i) a}- This effective en- 
tanglement presents a natural test for a quantum chan- 
nel, which outputs states p° ut for each input \ipi) from 
the test state ensemble. If we can verify entanglement 
between system A and the output states, the quantum 
channel did not act as a measure and re-prepare channel. 

A. Review of the expectation value matrix method 

The verification of entanglement between Alice and 
Bob bears two major difficulties. Firstly, the dimensions 
of H.a and TLb are, in general, different. The dimen- 
sion of system A is at most equal to the number N of 
pure test states. The dimension of the output system B 
may be much higher. The second difficulty arises from 
partial information. In most situations, it is challenging 
to perform state tomography on the output states, espe- 
cially in optical implementations (see, e.g., [HI)- ^ the 
measurements performed by Alice and Bob only supply 



partial knowledge, no unique state can be assigned to a 
set of measurement data. Instead, there will be a set of 
compatible states (an equivalence class). The task of ver- 
ifying entanglement is now mapped to checking whether 
this equivalence class contains a separable state. 

The (EVM) [23[ is a tool to perform exactly this task 
and it works for arbitrary dimensions and for a broad 
class of observables. Its construction depends on the 
test state ensemble and on Bob's measurement opera- 
tors. That is, the entries of the EVM are defined, as a 
function of the effective bipartite state pab, by 

[x(pA B )}ijki = Tr(p AB iJife <8> B]Bi). (1) 

The operator sets {Ai} and {Bj} must be chosen accord- 
ing to the measurement operators given by the protocol. 
For a binary input ensemble, the most compact set {Ai} 
is {|</>)(0|, \cf>) (1|}, where is a generic qubit state, which 
does not actually come into play due to the structure in 
Eq. (TT]). For each possible set of measurement outcomes, 
the EVM is designed to be a compact representation of 
the corresponding equivalence class of states. The mem- 
bership of separable states in this equivalence class is 
then tested by using the EVM in conjunction with a pos- 
itive, but not completely pos itive map [H|, such as the 
partial transposition map [26j | . By construction [22| , 

x{pab)>0 V pab- (2) 

Therefore, violation of the condition 

X(P T A A B ) > (3) 

is a sufficient condition for entanglement in pab- To 
check this condition, it must be possible to relate x(Pab) 
to the measurement outcomes stored in x(pab)- In the 
above choice of operators for Alice, x(Pab) * s related to 
x(pab) by a simple block transposition. 

If not all entries of x(pab) are experimentally acces- 
sible, the condition ([3]) can be checked with the help of 
semidefinite programming (27j . 

B. Pure squeezed test states 

One application of the EVM method was presented in 
Refs. [HJ and [H|], where the two coherent states |a) and 
|— a) are used to test a channel for quantum cryptography 
against intercept-resend attacks. Here, we investigate the 
similar task of testing a quantum memory with squeezed 
states, for which the same method can be employed. A 
convenient choice of test-states are the squeezed vacuum 
and the phase-rotated squeezed vacuum, 

\% n ) = S(r m )\0) (4) 
\i>?) = S(-r in )\0). (5) 

If two orthogonal quadratures x and p of the output light 
are measured, and the first and second moments of those 
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operators are recorded, an EVM can be constructed anal- 
ogously to Ref. [13]. 

It is at this point instructive to find a parameteriza- 
tion of typical measurement outcomes. We then show 
for which parameter regimes entanglement verification is 
possible. Since the operators x,p,x 2 , and p 2 are mea- 
sured on two possible output states, the measurement 
outcomes are characterized by eight parameters. How- 
ever, if the quantum device at hand does not introduce 
displacements, all first moments will be zero. Moreover, 
typical quantum channels are oblivious of the phase space 
orientation of the input states, so the output variances 
will be related by 



V&T p out (x) = Vax p o«t (p) 
Vai p out(p) = Var p o„ t (x). 

The EVM then takes the simple form 



/Var (i) b 1 + I ci c 2 + -§s\ 

bi - \ Var (p) c 2 - \s c 3 

c\ c* 2 + ±s Var (p) b 2 + | 

\c* 2 -y c* 3 6 2 -| Var (i)/ 



(6) 
(7) 



(8) 



where s denotes the overlap 
The operators Bj are chosen from the set {x,p}, since 
the first moments need not be included. The parameters 
bi e R and a € (D are not experimentally accessible, 
so that the condition §3§ must be checked by numerical 
evaluation of the resulting scmidefinite program. 

Now, two parameters, say Var p o U t(i) and Var p <pt(j5), 
suffice to specify a distinct set of measurement outcomes. 
Figure [1] shows the area of physically allowed measure- 
ment data in this parameterization. The gray area is ex- 
cluded by the Heisenberg uncertainty relation. Using the 
EVM method, we can identify parameter pairs which cor- 
respond to effective entangled states. These parameter 
pairs form the quantum domain, in which the quantum 
channel indeed operates in a non-classical manner. The 
boundary of the quantum domain shifts as the test-state- 
overlap s is varied. 

The plots show that the observation of squeezing 
alone, i.e., Var p o Ut (p) < 1/2, does not suffice. For higher 
degrees of squeezing in the test states, one must also 
observe more squeezing in the output states to verify 
operation in the quantum domain. 



C. Optimal measure and re-prepare channels 

The remainder of this section concerns those parame- 
ter pairs outside the quantum domain. Since the EVM 
method and indeed partial transposition form in general 
sufficient criteria for entanglement, but not necessary 
ones, points outside the quantum domain can fall into 
two categories: either they stem from effective entangled 
states which were not detected, or they are compatible 




FIG. 1: (Color online) Possible measured variances from two 
squeezed test states. For higher degrees of input squeezing, 
higher noise reductions must be observed to verify effective 
entanglement. The shaded area is excluded by the uncertainty 
principle. 



with a classical channel, i.e., a measure and re-prepare 
strategy. In the above parameter space, the latter cat- 
egory leads to a classical domain. By maximizing this 
classical domain, we can assess the strength of our test 
procedure. 

Mathematically, a measure and re-prepare str ateg y is 
equivalent to an entanglement-breaking channel 17], 



(9) 



where the operators iti form a positive operator valued 
measure (POVM), i.e., they are positive semidefinite and 
sum to the identity. Maximizing the classical domain 
now corresponds to the optimization of a suitable fig- 
ure of merit by varying both the measurement operators 
{iri} and the re-prepared states For the afore- 

mentioned case of coherent test states |a) and | — a), the 
figure of merit is the excess noise, or variance broadening, 
measured by Bob and the optimal entanglement-breaking 
channel minimizes this excess noise. This optimal chan- 
nel was found in Ref. p2| to be comprised of a minimum 
error discrimination POVM [28| and the resending of dis- 
placed squeezed states. For the case of two squeezed test 
states, the optimal entanglement-breaking channel is the 
one which minimizes Var p out (p) for each fixed value of 
Var pr (i) (c.f. Fig.©. 

Let us first consider the re-prepared states {(V'i)}- 
Since the figure of merit is a variance, the re-prepared 
states must be minimum-uncertainty states. Further- 
more, since Bob expects all first moments to be zero, 
displacements in {iV'i)} have the effect of increasing the 
quadrature variances (see, e.g., [29(). This leaves the 
class of squeezed vacuum states, with the squeezing axis 
aligned with the axis of Bob's quadrature detection. 
Therefore, each re-prepared state can be parame- 
terized by a single real parameter, namely, its squeezing 
parameter r;. Let us now turn to the POVM {iii]. Since 
there are only two pure test states, each POVM element 
is described by a two-dimensional matrix. Assuming that 
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Var (x) 

FIG. 2: (Color online) Optimal measure and re-prepare strat- 
egy for two squeezed test states. The boundaries of the clas- 
sical and the quantum domains coincide. Optimal measure- 
ments in region (A): minimum error discrimination; region 
(B): unambiguous state discrimination 



while Var p crut (p) remains finite, explains the straight line 
with slope (approaching) zero. Such a strategy involves 
the unambiguous discrimination of the two input states 
[3fj| . When the measurement outcome is inconclusive, 
vacuum states are re-prepared, and whenever is 
identified unambiguously, an infinitely squeezed state is 
re-prepared. This way, Var p out (x) tends to infinity while 
Var p out (p) remains finite. 

With the combination of minimum error measurement 
and unambiguous state discrimination, we can show an- 
alytically that the EVM method witnesses all entangled 
states which are detectable by the given measurements. 
The resulting quantum domain cannot be enlarged by 
any other verification method. This is in accordance with 
the results for two coherent test states |a) and | — a). 



the number of POVM elements is small, the optimal so- 
lution can easily be found numerically. Explicitly, the 
figure of merit is given by 

Var, r (p) = i^tr(7r l p™)exp(-2r J ), (10) 



which we wish to minimize for each fixed value of 
Var p out (x) under the constraints 



7T, >0, 

5>< = i. 



(11) 

(12) 



The results of a numerical optimization are shown in 
Fig. [2j with the number of POVM elements set to 
four. We see that this optimized measure and re-prepare 
strategy reaches the boundary of the quantum domain. 
Therefore, any experimental outcome in our parameteri- 
zation is unambiguously sorted into one of three domains: 
the unphysical domain, the quantum domain as deter- 
mined by the EVM method, and the classical domain. 

It is natural to ask whether this optimal POVM cor- 
responds to minimum error discrimination, in analogy 
to the coherent-state-protocol. It turns out that this 
is true only partially. The boundary of the quantum 
domain shown in Fig. [2] is a concatenation of a curved 
part (A) and a straight line with slope zero (B). This 
behavior arises from the fact that the first eigenvalue of 
the EVM to become negative is different in regions (A) 
and (B). Therefore, the optimal measure and re-prepare 
strategies and in particular the POVMs are also likely to 
be different in both regions. It turns out that the optimal 
measurement for region (A) is indeed the minimum 
error discrimination. For the measurement in region 
(B), we recall that the measured variances are linear 
functions of Pq"*, so that any convex combination of 
two measure and re-prepare strategies results in another 
valid entanglement-breaking channel. Therefore, finding 
one particular strategy for which Var p ™t (£) diverges, 



III. MIXED STATE SOURCES 

In this section, we consider sources of mixed test states. 
In a communication context, this mixing operation may 
be useful and simplify the transmission of secret messages 
[3l| . More common, however, are situations where it may 
be desirable to use pure states, but not experimentally 
feasible. 

Any source of mixed quantum states p™ 1 can be 
thought of as a source of pure states |V>| n ); followed by a 
mixing operation. This picture shows a natural way to 
include mixed-state sources in the EVM-method. For the 
entanglement-based description and the construction of 
the EVM, we can consider the pure states \ipl n ), and at- 
tribute the mixing process to the quantum channel to be 
tested. In this way, the method from the previous section 
can be used. However, the mixing process will introduce 
an additional degradation of the channel. Therefore, any 
quantum benchmark derived in this way will be valid, 
but not tight, since it actually is a benchmark on the 
concatenation of the mixing process and the action of 
the quantum channel. 



A. Mixed test states and purifications 

What is a better way to describe the mixed-state 
sources? Since the EVM-method works very well for 
sources of pure states, we would like to build on the lan- 
guage developed in the previous section. We start by 
fixing some notations. Let the test state ensemble be 
given by the probabilities pi and the corresponding test 
states p\ n . We can decompose each test state into pure 

states as pf = £\ qf\j^)&% 

Considering each mixed state emitted from the source 
as a statistical mixture of pure states, we see that 
the source actually emits pure states from the set 

,(<h Fig.|3(a)1. 



w ith probabilities {pi qf 1 } 



An 

entanglement-based description of the source is now 
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(a)mixture of pure states (b)part of a purification 

FIG. 3: Each mixed state emitted from the source can be 
regarded as a mixture of pure states, or as part of a higher- 
dimensional pure state (purification) 

straight-forward. Problems arise when it comes to the 
construction of the EVM, since Bob's measurements can- 
not resolve which pure state in the decomposition of a 
test state was sent. In other words, the measurement 
outcomes may only be conditioned on the index i, but 
not on the index j. This makes the direct construction 
of an EVM impossible. 

There is, however, an alternative way to connect mixed 
states and pure states. Every mixed state p can be 
thought of as a part of a higher-dimensional pure state 
I?/;), which is called a purification of p [321 ] . Hence, a 
possible simplified description of the actual mixed-state 
source is a source which emits purifications of the test 
states. This is illustrated in Fig. 3(b) The test-state 



source, which emits mixed states on system B is ap- 
proximated by a source which emits purifications \tpi) bc, 
where the dimension of system C is at most equal to the 
dimension of system B. For simplicity, we will again 
concentrate on the smallest non-trivial input ensemble, 
which contains only two test states, each occurring with 
probability 1/2. In an entanglement-based description, 
the resulting effective entangled state is 



1 



\r rc ) = -j=(\0)Mbc + \1)a\A)bc), 



(13) 



where the \ipijBC & r e purifications of the actual test 
states pf 1 , i € {0,1}. 

An EVM can now be constructed in direct analogy 
to the pure-state-case of the previous section. The fact 
that Bob does not have access to the auxiliary system 
is easily incorporated by appending lc to each operator 
which enters the EVM. In other words, in the general 
form of a three-party EVM 



\x(pABc)]ijklmn = ^{pABC A\Ai ( 



>b]b„ 



<C n ), (14) 



the set {Cfc} has only one member, namely lc- This 
way, the auxiliary system and the use of purifications do 
not increase the size of the EVM when we compare it to 
the pure-state-case. The dependence of the EVM on the 
test-state-overlap s is changed to a dependence on the 
overlap of the two purifications \4>o)bc and \4>i)bc- 



B. Optimal purifications 

It should be noted that a purification is not unique, 
i.e., there is a freedom in the choice of \ipi)BC- It turns 




FIG. 4: Source-replacement scheme for a source of mixed 
qubit states. 



out that the entanglement criterion works best if \iJjo}bc 
and bc are chosen such that their overlap takes its 
maximum value. From Uhlmann's theorem [331 ], we know 
that 

max IW«|=W,/-r), (15) 
lw,|V , i) 

where the maximization runs over all purifications 
of the states p\ n , and F denotes the fidelity 



F{pa,pi) = trW ^a/poPiVPo". 



(16) 



As mentioned above, a source which emits the purifi- 
cations \tjjijBC is an approximation to the true test-state 
source. This approximation leads to weakened bench- 
marks, since the two purifications can be distinguished 
more easily than the mixed test states, which facilitates 
measure and re-prepare strategies. Ideally, only system 
B is emitted and the purifying system C is retained in 
the source [see Fig. |3(b)] . In other words, Alice has full 
access to the reduced density matrix Pac- However, the 
above construction of the EVM does not contain the full 
matrix pac, but only pa- 



C. Further improvements 

In order to witness all of the detectable entanglement, 
the information stored in the purifying system must be 
included in x(pabc)- A natural way to achieve this is 
through measurements on system C. The correspond- 
ing measurement operators {Cfc} can then be included in 
the EVM, c.f. Eq. (fT4"]) . If a tomographically complete set 
{Cfc} is chosen, pac will be fully contained in x(pabc)- 
However, care must be taken with this approach. The 
test states must not be conditioned on measurement out- 
comes of the system C. This would lead to a different 
test state ensemble. Also, expectation values on system 
C are only available after tracing over Bob's system, since 
they reflect the knowledge of the reduced density matrix 

PAC- 

Below, we describe three physical setups: One qubit 
protocol, one setup with squeezed thermal states as test 
states and one setup with displaced thermal states. 

Qubit states. We begin by considering an example with 
qubit test states, which will serve to set the notations 
and to demonstrate the use of purifications. The setup 
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FIG. 5: (Color online) Results of the entanglement verifi- 
cation for qubit states. The verifiable quantum domain for 
mixed states shrinks (dashed line) unless the optimal purifi- 
cations are used (circles) 

is shown in Fig. Q] Inside the test-state source, an entan- 
gled state \4> stc )ab is created, such that two pure qubit 
states |V'o rl )s an d IV'i ra )-B can be prepared by projective 
measurements on system A. It is convenient to choose 
them 

|<> s = cos(~)|0>B+sm(~)|l) B , (17) 

\^) B = C o S (^)\Q) B -sm( 9 -)\l) B . (18) 

These states are mixed by the interaction with an aux- 
iliary qubit system C through a conditional phase gate. 
This auxiliary system is initially in the superposition 
(|0) + |l))/v / 2- By adjusting the phase (j), the mixed- 
ness of the conditional states in system B can be varied 
continuously. 

In this setting, = corresponds to pure test states. 
For this case, no purifying system is required. The re- 
sults of the entanglement verification are shown in Fig. [5] 
Here, Bob can measure all three Pauli operators and his 
observations are modeled by a depolarizing channel, 

P ^=KV(^)=P^ + {1-P)P?- (19) 

The figure shows the maximum value of the depolariz- 
ing parameter p for which effective entanglement is still 
detected, for different values of the angle 8. We observe 
that more channel noise can be tolerated when the signal 
states overlap more, similar to the observations for two 
phase-conjugated coherent states. 

What happens when the test states become mixed, i.e., 
the conditional phase gate is in place with 0=^0? The 
joint states \ipi)BC of systems B and C after this gate 
will not be of product form anymore, and they are pu- 
rifications of the mixed test states. The dashed line in 
Fig. O shows that the verifiable quantum domain shrinks. 
However, the unitary freedom in the choice of the purifi- 
cations has not been exploited yet. If we choose those 
two purifications with the maximum overlap, more effec- 
tive entangled states are detected (circles in Fig. [SJ. In 



fact, the detected quantum domain is the same as for 
pure test states. 

Finally, we can include the information contained in 
the purifying system C, by performing measurements on 
it. The Pauli operators are an obvious choice of measure- 
ment operators, since, together with the identity, they 
form an operator basis for Tic- I n practice, this amounts 
to including the set {Ck} = {lei fi, &y} in the EVM. It 
turns out that this larger EVM cannot detect any addi- 
tional states, or in other words, the approximation with 
purifications is an equally strong criterion. 

This a particularity of this example with qubits, which 
will be examined in Sec. IIVI The next two examples 
show that the inclusion of the information stored in pac 
can indeed improve the EVM-method. 

Squeezed thermal states. Squeezed light is of funda- 
mental interest due to its non-classical character, which 
shows for example in sub-Poissonian photon number dis- 
tributions and in uncertainties below the level of vacuum 
fluctuations. The generation of quadrature squeezed light 
by processes such as optical parametric amplification has 
become a standard technique. The use of two squeezed 
vacuum states to probe a quantum channel was described 
in the previous section. An extension of the method to 
mixed states is essential, since squeezed states become 
thermalized under the influence of losses. The resulting 
test states are 

^"-^-tEI^XT S(r m )\n)(n\S\r m ) (20) 
n + 1 \ n + 1 J 

n—0 v / 
oo / _ \ n 

^ l = ^TTEMrT SHr m )\n)(n\S(r m ), (21) 
n + 1 \n + 1 / 

n— 

i.e., they are characterized by their mean photon num- 
ber n and the squeezing parameter r™ (see, e.g., (29|). 
The source replacement and the entanglement verifica- 
tion work very similarly to the above example with qubit 
states. The main differences are the following: While the 
qubit test states act on a two-dimensional Hilbert space, 
independently of the degree of mixing, two squeezed ther- 
mal states act on an infinite-dimensional Hilbert space, 
and each test state has full support. Furthermore, full to- 
mography on the output states would result in an infinite- 
dimensional EVM, so, in analogy with the previous sec- 
tion, we consider detection of the first and second mo- 
ments of the quadratures x and p. For the typical mea- 
surement outcomes mentioned above, i.e., both quadra- 
ture operators have zero expectation values and the vari- 
ances are related through Eqs. ^ and (O, the EVM 
for mixed states is of the same form as for pure states 
[Eq. JSJ)]. Here, the operator set {Ck} in the general 
form (|14[) has only one member, namely, 1q- 

Results of the entanglement verification are shown in 
Fig. [6] These results strongly resemble those from the 
previous section, where pure squeezed test states were 
considered. Indeed, the only difference is that the over- 
lap of the pure test states is replaced by the fidelity of 
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FIG. 6: (Color online) Entanglement verification from dif- 
ferent EVMs. Solid line: no information on the purifying 
system; dashed line: one generalized spin operator on the pu- 
rifying system; dotted line: three generalized spin operators; 
dashed-dotted line: five generalized spin operators 



the mixed test states. We therefore note that the EVM- 
method can be extended to include sources of mixed 
states without introducing additional complexity. But, 
as mentioned above, we are merely approximating the 
source of two mixed test states by a source of two purifi- 
cations of the test states. 

We now proceed to include more of the information 
stored in pac m the EVM. This is done by enlarging 
the set of operators {Ck} which enter the construction 
of x(pabc)- For the numerical implementation, it is 
convenient to choose unitary operators, and ideally, the 
set becomes large enough to form an operator basis 
for the space of operators on Tic- However, since the 
dimension of this space is infinite, we are forced to 
restrict ourselves to a small subset. It is also unknown 
which operators to choose in order to achieve the biggest 
improvement of the entanglement criterion. We therefore 
choose the operators from a generic set, namely the 
generalized spin operators (see, e.g., [Hj]), which can 
be regarded as a generalization of the Pauli operators 
to higher dimensions. Each additional member in the 
set {Ck} enlarges the dimension of the EVM, and we 
have evaluated the entanglement criterion with up to 
five additional unitary operators on system C. Here, 
we worked in a truncated Fock basis with maximum 
photon number n max = 16. The corresponding results 
are shown in Fig. [5] It can be seen that the addition 
of operators on the purifying system C does indeed 
improve the criterion. The more operators are included, 
the more entangled states are detected, but the size of 
X grows and with it the number of free parameters. 
Therefore, there is a trade-off between the strength of 
the entanglement criterion and its complexity. It will, 
however, become clear in the next section that the 
curves shown in Fig. [6] are already close to the border of 
the classical domain. 

Displaced thermal states. Finally, we briefly examine 
the case of two displaced thermal states as test states. 
Experimentally, coherent states can be generated essen- 
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FIG. 7: (Color online) Quantum domain for displaced ther- 
mal test states. (a) the quantum domain increases in 
size for mixed states, (b) improvements of the entangle- 



ment verification by adding two unitary operators to the 
EVM (dashed) and optimized measure and re-prepare strat- 
egy (dashed-dotted) 



tially noiselessly. Still, the addition of noise to the test 
states may be of advantage to the protocol at hand. 
Moreover, this second example in the continuous- variable 
domain will supplement the results for squeezed thermal 
states and allow us to draw a comparison. 

In terms of Fock states, the two test states are given 

by 



1 



Po 



Pi 



n=0 



tE(^Ti) D(a)\n){n\D^a) (22) 



n=0 



tE(^Ti) D(-a)\n){n\DH-a). (23) 



Here, D(a) denotes the displacement operator [29( and 
we take a to be real. These two test states pass through 
the quantum channel and, as before, homodyne detec- 
tion is performed on the output states. We will, however, 
choose a different parameterization of the possible mea- 
surement outcomes in order to be consistent with [22l. l23l] . 
Specifically, the observed variances of the x and p quadra- 
tures are typically equal, and losses act as a simple down- 
scaling of the observed first moments. The results for 
pure test states [H, [23[ are displayed in Fig. 7(a) 
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As before, a good approximation to the source of mixed 
test states is a source of corresponding purifications. Re- 
sults of the entanglement verific ation using the optimal 
purifications are shown in Fig. 7(a) Remarkably, the 



verifiable quantum domain is enlarged compared to the 
setting with pure states. To be precise, two coherent test 
states \a) and | — a) lead to a smaller quantum domain 
than the same coherent states with added noise. 

The strength of the entanglement criterion can be im- 
proved by including information on the purifying sys- 
tem C . As for the above example with squeezed thermal 
states, we enlarge x(pabc) with unitary operators Uc 
from the set of generalized spin operators. The dashed 
line in Fig. |7(b)| shows the results for two added uni- 
tary operators. Evidently, the quantum domain is sig- 
nificantly enlarged. But how close are these results to 
the boundary of the classical region, i.e., to those mea- 
surement outcomes compatible with a measure and re- 
prepare channel? The results of such a channel arc indi- 
cated by the dashed-dotted line in Fig. |7(b)| Here, the re- 
prepared states are displaced squeezed vacuum state, and 
the measurements are optimized numerically. A more 
detailed discussion of optimized measure and re-prepare 
channels is given in the next section. 

We note that the results for displaced thermal states 
and squeezed thermal states are qualitatively very simi- 
lar. The approximation of the true test-state source by 
a source of purifications leads to simple and reasonably 
strong entanglement criteria. In both cases, the addi- 
tion of a small number of unitary operators which act on 
system C pushes the boundary of the detected quantum 
domain close to the boundary of the classical domain. 



IV. OPTIMIZED 
ENTANGLEMENT-BREAKING CHANNELS 

As before for the pure test states, it is of interest to 
find the border of the quantum domain the converse way, 
i.e., by optimizing corresponding entanglement- breaking 
channels. According to Eq. ©, an entanglement- 
breaking channel is characterized by the POVM {71^} and 
the re-prepared states {IV'i)}- In this section, we will de- 
rive optimized measure and re-prepare strategies for the 
qubit setup and the squeezed state setup described above. 

For pure qubit test states, it is again the minimum er- 
ror POVM which leads to the optimal strategy. If the 
re- prepared states are chosen of the same form as 

the original test states, the noise added by the measure- 
ment errors can be expressed as the action of a depo- 
larizing channel. This allows a comparison between the 
boundaries of the quantum domain as shown in Fig. [5] 
and the classical domain, and we indeed find that both 
boundaries coincide. Therefore, the EVM method pro- 
vides a sharp boundary of the quantum domain for pure 
qubit test states, as one might expect from the results of 

Sec. nu 

We recall that the mixing process inside the test-state 



source did not shift the boundary of the quantum domain 
for qubit states. In fact, this boundary can be reached 
by an optimized entanglement-breaking channel for any 
degree of mixing of the test states. The optimal POVM 
is that of minimum error discrimination. This may seem 
surprising, because the error probability, which depends 
on the fidelity (or the overlap) of the test states when 
they are pure, is a function of the trace distance between 
the test states when they are mixed. As it turns out, the 
well known relation between fidelity and trace distance 



D tr (p,a) < y/1 -F(p,a) 



(24) 



is fulfilled with equality if p and a are qubit states with 
equal degree of mixing. This can be shown by direct 
calculation with the help of the results in Ref. [35j . 
Hence, for the investigated qubit system, mixed signal 
sources can be incorporated in the EVM-criterion with- 
out adding complexity to the problem. The resulting 
boundary of the quantum domain is tight, i.e., the 
quantum domain determined by the EVM-method 
reaches the boundary of the classical domain. 

Squeezed thermal states. In Sec. HH we derived the 
optimal measure and re-prepare strategy for an ensem- 
ble of pure test states. In part, we can use these re- 
sults when considering mixed test states. For instance, 
it was argued that the re-prepared states {JV'i)} should 
be drawn from the set of squeezed vacuum states. Natu- 
rally, this is also true for mixed test states. The optimal 
POVM, however, is of a very different structure. While 
two pure test states span a two-dimensional space and the 
POVM elements are consequently represented by two- 
dimensional matrices, the squeezed thermal states act on 
an infinite-dimensional Hilbert space. Nevertheless, the 
optimal POVM can be found analytically. 

We recall the above figure of merit for squeezed input 
states: The optimal entanglement-breaking channel min- 
imizes Varo (p) for a fixed value of Varo (x) . It is instruc- 
tive to start with the simplified problem of finding the 
global minimum of Var p o«t (p) while ignoring Var p <^t (£) 
for the moment. The only constraints on the minimiza- 
tion are that p° ut are the results of an entanglement- 
breaking channel acting on the input states p™, and 
Eqs. (|6]) and ([7]) must be fulfilled. In fact these last 
two constraints are automatically satisfied if we consider 
entanglement-breaking strategies which obey a certain 
symmetry. If {n.i, \ipi)} is a particular strategy which 
fulfills Eqs. (O and ([7]), then the phase-rotated strategy 
{UziTiUl , leads to the same observations. Here, 



Ujl = exp(-i-n) 



(25) 



is a phase shift operator. An equal mixture of both 
strategies will then automatically satisfy the constraints 
© and HI). We find 
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\ i 2 



(26) 



For each term in the summation, we can find the opti- 
mal degree of squeezing, r^, and insert it into the above 
expression. We find 

Var pr (P) = \Y, JtT(n iP ™)tT(n iP ?). (27) 

i 

For the remaining minimization over the POVM ele- 
ments, we borrow a result from Ref. [36| . namely, 



^trfa^trfari") = F(j$,fr), (28) 



where F denotes the Uhlmann fidelity defined in Eq. [THl 
In summary, the smallest variance achievable by a mea- 
sure and re-prepare strategy is given by half the fidelity 
of the two test states (see Fig. [6]) . This is in agreement 
with the results of the entanglement verification. The 
boundary of the quantum domain is therefore tight at 
this minimum point. But is it tight elsewhere, i.e., for 
different values of Var p o U * (£)? 

To answer this question, we introduce the additional 
constraint 



Vaxpaut (x) 



(29) 

to the optimization problem, and a plot like Fig. [8] is 
obtained by varying c in steps and minimizing Var p ou* (p) 
at each step. This problem can be solved using Lagrange 
multipliers, i.e., we can minimize the Lagrange function 

C = Var pr (p) + A(Var pr (4) - c), (30) 

which is now unconstrained, but the minimization runs 
over {7Ti}, {ri\, and A. Since both variances are linear in 
Pq™ and in p™, we can write the above expression as 

£= l(^ (E " 2ntr( ^ [p ° l + M " ]) 



W r M^[pT + ^\)) - Ac. 



Defining the new states p = (p™ + Apj n )/(1 + A) and 
pi = (p" 1 + A/9q 1 )/(1 + A), and minimizing over the Vi, 
we arrive at 

C = 1 ~t ^ y/tr(7rjpo)tr(7Tipi) - Ac (31) 



i£ = 



2 

1 + A 



F (po,Pi) - 



(32) 
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FIG. 8: (Color online) Optimal measure and re-prepare strat- 
egy for squeezed thermal test states. Solid: entanglement ver- 
ification with no information on the purifying system; dashed- 
dotted: five additional unitary operators on system C in the 
EVM; dashed: optimized measure and re-prepare strategy. 



In the last step, the relation (|28|) was used again. Now 
we can minimize this function with respect to A, subject 
to the original constraint (|29|) . Although this must still 
be done numerically because of the fidelity function, the 
problem is considerably simpler than a complete opti- 
mization over the POVM elements. 

As it turns out, this method does not work for all val- 
ues of c. Drawing the analogy to the pure-state case and 
Fig. [U only the optimal strategy in region (A) can be 
explained thus. In region (B), the Lagrange parameter 
A becomes negative, and the states po and p\ turn un- 
physical. For pure test states, the optimal measurement 
in region (B) depended on unambiguous state discrimi- 
nation. Such a measurement strategy is not directly ap- 
plicable [13] to two squeezed thermal states, since they 
have equal support. It is still possible to construct a mea- 
sure and re-prepare strategy which reaches the straight 
line bounding the quantum domain. A specific strategy 
which approximates the optimum arbitrarily well is pre- 
sented in Appendix [AJ 

We now see that the entanglement verification method 
does not detect all states in the quantum domain, but 
increasing the size of x(pabc) and including more infor- 
mation about system C brings a notable improvement. 
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CONCLUSION 



In conclusion, a generic method to test quantum chan- 
nels was presented. The method requires very few ex- 
perimental resources and is therefore directly applicable 
to quantum memory or teleportation experiments. To 
accommodate tests with mixed states, we extended the 
concept of effective entanglement. Verification of such 
effective entanglement was investigated using the expec- 
tation value matrix in conjunction with partial transpo- 
sition. A strong, yet simple extension of the expectation 
value matrix method for pure-state sources was proposed, 
where a source of mixed test states is approximated by a 
source of purifications of these test states. This approxi- 
mation could detect all entanglement for the considered 
qubit state example. For two continuous-variable set- 
tings, where displaced thermal states and squeezed ther- 
mal states were considered, using purifications of the true 
test states lead to strong, but sub-optimal criteria. A 
procedure to improve the EVM method for these cases 
was proposed and conjectured to lead to optimal entan- 
glement criteria. Furthermore, optimized measure and 
re-prepare strategies were considered as the counterparts 
of entanglement-preserving channels, and as a means of 
probing the strengths of the entanglement criteria. 

The results obtained for modeled channels suggest that 
the derived benchmarks are reachable by current experi- 
ments. 
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APPENDIX A: USD-TYPE STRATEGY FOR 
MIXED STATES 



Here, we derive a measure and re-prepare strategy for 
two squeezed thermal states as test states, for the region 
(B) (Fig. [9]). Again, we use the fact that the measured 
variances are linear functions of pg u * , so that any convex 
combination of two measure and re-prepare strategies re- 
sults in another valid entanglement-breaking channel. It 
therefore suffices to consider the slope of the straight line 
which joins the point p in Fig. [5] and a the point aris- 
ing from an entanglement-breaking channel in the limit 
Var p o Ut (i) -> oo. 

Let the coordinates of p be denoted by (v x ,v p ), then 
the slope of the straight line joining p and any other point 



P 



region (A) 



region (B) 



0.7 0.8 0.9 

Var (x) 



FIG. 9: Schematic plot of the two different optimal measure 
and re-prepare strategies. 



is given by 



Var pr t (p) - v p 
Var p o„ t (x) - v x ' 



(Al) 



where Var pOU t (p) is given by Eq. (|26p and the variance of 
x is defined accordingly. 

The strategy which minimizes this slope is essen- 
tially an approximation to unambiguous state discrimi- 
nation. We define a three-element POVM with IIo, IT = 
L r Jy 2 rioL r - 7r /2 , and ft?, and corresponding re-prepared 

states l^o) = \ r ), = | — t), and \ip-?) — |0). Con- 
sequently, 



Var prt (p) = i(tr(p»7T )e- 2 '- + tr(pl"7r )e 2 '' + tr(p« n 7r ? )) 



1 



Var pr (£) = -(tr(p-7r )e^ + tr^TTo)^ + tr(p™7r ? )). 



Inserting into Eq. (lAlj) and taking the limit r — > oo, we 
find 



lim m 



tr(pf7ro) 



(A2) 



All that is left to do now is to find a POVM element ttq 
for which the slope to tends to zero. Setting ttq = \a)(a\, 
where \a) is a coherent state, achieves this. In this case, 
to becomes a quotient of Husimi Q- functions: 



\<x\Pi\<x) 
(a\p \a) 



Qo{a) 



(A3) 



The Q-function Qi(a) of the squeezed thermal state pi, 
i G {0, 1}, is given by [38j 



1 



7r v / (l + ^) 2 -|S l | 2 



eM(A-l)\a\ 2 -l(B ia * 2 +B*a 
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with the definitions Inserting these expressions into Eq. (|A3|) , we find 



A = n+(2n + l)smh 2 {r)=A 1 (A4) expl-^a* 2 + a 2 )} - . . 

B = -(2n + 1) sinh(r) cosh(r) = -B 1 (A5) m = ^^{a* 2 + a 2 )] = e ^~ B ^ a * 

n(n + 1) (A6) 
n 2 + (n + |)(1 + cosh(2r)) For increasing a € R, this final expression tends to zero 

/- . i\ s j n k{2 r ') exponentially fast. Therefore, it is possible to reach the 

B = — — -; \— — = —B\. (A7) boundary of the quantum domain in Fig. [3] for the whole 

n 2 + (n + + cosh(2r)) parameter domain. 
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